function R = lbppe_features(image)
R = [processOne(image, 1), ...
    processOne(image, 2), ...
    processOne(image, 3), ...
    processOne(image, 4), ...
    processOne(image, 5)];
end
   
function f = processOne(image, k)
r = compute_residual_image_succ(image, k);
f = [lpq(r, 3, 0, 1, 'nh'), ...
    lpq(r, 5, 0, 1, 'nh'), ...
    lpq(r, 7, 0, 1, 'nh')];
end
